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We address the issue why the phase diagrams for quasi-one-dimensional systems are rather simple, 
while the renormalization group equations behind the scene are non-linear and messy looking. The 
puzzle is answered in two steps - we first demonstrate that the complicated coupled flow equations 
are simply described by a potential V(hi), in an appropriate basis for the interaction couplings 
hi. The renormalization-group potential is explicitly constructed by introducing the Majorana 
fermion representation. The existence of the potential prevents chaotic behaviors and other exotic 
possibilities such as limit cycles. Once the potential is obtained, the ultimate fate of the flows are 
described by a special set of fixed-ray solutions and the phase diagram is determined by Abelian 
bosonization. Extension to strong coupling regime and comparison with the Zamolodchikov c- 
theorem are discussed at the end. 



I. INTRODUCTION 

Quasi-one-dimensional (Q1D) systems have attracted 
extensive attentions from both experimental and theo- 
retical aspects^. Due to the low dimensionality, strong 
quantum fluctuations often give rise to surprising behav- 
iors, which are rather different from our intuitions built 
in higher dimensions, ft is then exciting to explore var- 
ious exotic phenomena, such as spin-charge separation 2 , 
unconventional electron pairing^ and so on, in the cor- 
related Q1D systema 5 i 6 i 7 i 8 iV . Note that, in addition to 
the ladder compounds, carbon nanotubes, nanoribbons 
and quantum wires, after integrating out fluctuations at 
higher energy, are also described by the Q1D theory. 
Therefore, not only posting a challenging task for aca- 
demic curiosity, the understanding of the Q1D systems 
now becomes crucial as the technology advances to the 
extremely small nanometer scale. 

Although our understanding of the strictly one- 
dimensional systems is greatly benefitted from exact solu- 
tions, it is often found that the Q1D systems are not sol- 
uble analytically. Furthermore, since the hopping along 
the transverse direction is relevant, the physics phenom- 
ena for Q1D systems can be dramatically different from 
the strictly ID systems, such as electrons with mutual 
repulsive interactions are found to form unconventional 
"d-wave" Cooper pairs in the ladder systems. Therefore, 
the most reliable approach to clarify the competitions 
among various ground states is the renomalization group 
(RG) analysis. Since the number of allowed interactions 
are large, the derivation of the flow equations for all cou- 



plings under RG transformations becomes formidable. 
On the other hand, both numerical and analytical ap- 
proaches seem to indicate rather simple phase diagrams 
for generic QfD systems. 

The simplicity of the phase diagram can be partially 
understood by the widely used scaling Ansatz of the cou- 
plings in weak coupling2*2*2*ii, 
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where Gi are order one constants satisfying the non- linear 
algebraic constraints (discussed later) and Id is the diver- 
gent length scale where the flows enter the strong cou- 
pling regime. The Ansatz was motivated by the numeri- 
cal observation that the ratios of renormalized couplings 
reach constant, as long as the bare interaction strength 
is weak enough. However, it is still puzzling why the 
phase diagrams, generated by many coupled non- linear 
flow equations, do not reflect the same level of complex- 
ity. In fact, for a complicated system with many cou- 
plings, the coupled non-linear differential equations are 
likely to produce chaotic flows generically. Even if the 
flow is not chaotic, it might as well rest on limit cycles. 
This peculiar possibility for quantum systems was ad- 
dressed by Wilson and collaborators in a recent paper—. 

So the question remains: Why are the phase diagram 
and the RG flows so simple (without chaotic flow or limit 
cycle) in QID systems? We found the question can be 
answered by combining the weak-coupling RG analysis 
together with the non-perturbative Abelian bosonization 
technique. Note that the phases of the correlated ground 
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state often lie in the strong coupling, so weak-coupling 
RG alone can not pin down the phase diagram. It is the 
powerful combination of the perturbative RG analysis 
and the non-perturbative bosonization technique which 
deliver the desired answer here. 

Let us start with the weak-coupling RG. In the low- 
energy limit, the Q1D systems involve Nf flavors of in- 
teracting fermions with different Fermi velocities, where 
Nf is the number of conducting bands. Although con- 
strained by various symmetries, the number of allowed 
interactions is tremendously large as Nf increases. In 
general, the RG equations to the lowest order are al- 
ready very complicated - not to mention solving them 
analytically. 

However, quite to our surprises, we found that, at one- 
loop level, the RG flows cab be derived from a potential, 
i.e. the coupled non-linear flow equations can be cast into 
this elegant form by an appropriate choice of coupling 
basis, 

dh t _ dv(hj) 

dl dhi ' { ' 

where V(hj) is the RG potential. We emphasize that 
this is only possible after a unique transformation of the 
couplings, hi(l) — Lijgj(l) (up to a trivial overall fac- 
tor for all couplings), where is some constant matrix. 
The existence of the potential, which requires the coeffi- 
cients in the RG equations to satisfy special constraints, 
also provides a self-consistent check on the RG equations 
derived by other approaches. 

The flows of Eq. J2J) in the multi-dimensional coupling 
space can be viewed as the trajectories of a strongly 
overdamped particle moving in a conservative potential 
V(hi). Note that the change of the potential V(hi) along 
the trajectory is always decreasing, 

dv = d^dfn = _ f dhA 2 

dl dhi dl \dl J ~ ' U 

where summation over the index i is implicitly implied. 
Therefore, it is obvious that the function V(hi) never in- 
creases along the trajectory and is only stationary at the 
fixed points where dhi/dl = 0. Thus, the RG flows have 
a simple geometric interpretation as the trajectory of an 
overdamped particle searching for potential minimum in 
the multi-dimensional coupling space. 

This simple geometric picture rules out the possibili- 
ties of chaos and the exotic limit cycles in Q1D systems. 
The ultimate fate of the flows would either rest on the 
fixed points or follow along the "valleys/ridges" of the 
potential profil e - 13 ! 14 to strong coupling. Since there is 
only one trivial fixed point at one-loop order, most of the 
time, the flows run away from the non-interacting fixed 
point. Starting from weak enough bare couplings, the ul- 
timate fate of the flows is dictated by the asymptotes of 
the "valleys/ridges" of the potential profile. It provides 
the natural explanation why the ratios of the renormal- 
ized couplings reach constant in numerics. That is to say, 



the existence of RG potential implies that the ultimate 
fate of RG flows must take the scaling form described in 
Eq. ^ Detail properties of these asymptotes, referred as 
fixed rays, will be discussed in later section. 

Since the ultimate fate of RG flows is described by the 
simple Ansatz in Eq. ^ the specific ratios of couplings 
simplify the effective Hamiltonian a lot. Making use of 
the Abelian bosonization, one can determine which sector 
acquires a gap, triggered from the weak-coupling instabil- 
ity. The phase of the ground state is then determined by 
watching which fixed ray (asymptote) the flows go closer 
to. Because there are only limited solutions of the fixed 
rays, the phase diagram is rather simple. Therefore, by 
combining the powerful techniques of weak-coupling RG 
and Abelian bosonization, we pin down the reason be- 
hind the simple-looking phase diagram out of the messy 
non-linear flow equations. 

In fact, the combination of weak-coupling RG and 
Abelian bosonization goes beyond the usual mean-field 
analysis and is crucially important when there are more 
than one competing orders^. For instance, Lauchli, Hon- 
ercamp and Rice recently studied the so-called "d-Mott" 
phase in one and two dimensions, where antiferromag- 
netic, stagger-flux and d-wave pairing fluctuations com- 
pete with each other simultaneously. The conclusion 
drawn from the numerical density-matrix RG in strong 
coupling agrees rather well with predictions made from 
one-loop analysis in weak coupling. This lends support 
to the powerful combination of weak-coupling RG and 
Abelian bosonization approach for strongly correlated 
systems. 

Since the method of bosonizing the fixed rays is already 
developed in previous papersa^S, we would concentrate 
on the novel existence of RG potential in this paper. In 
particular, we would construct the RG potential explic- 
itly. The details of the bosonization and numerical re- 
sults will be deferred for future publication. The remain- 
ing part of the paper is organized in the following: In 
Sec. II, the criterion for potential flows and the notation 
of fixed rays are explained. In Sec. Ill, we write down 
the effective Hamiltonian for generic Q1D systems and 
elaborate all possible interaction vertices in weak cou- 
pling. In Sec. IV, the Hamiltonian is re-written in terms 
of Majorana fermions and the RG potential is explicitly 
constructed. Finally, we close the paper by the section 
on discussions and the summary of the main results. 



II. CRITERION FOR POTENTIAL FLOWS 

To prove the existence of the RG potential, it is help- 
ful to study the general feature of one-loop RG equations 
of the Q1D systems first. In weak coupling, the most 
relevant interactions are the marginal four-fermion inter- 
actions, described by a set of dimensionless couplings gi. 
The RG transformation to the one-loop order is described 
by a set of coupled non-linear first-order differential equa- 
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tions 
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where the coefficients Mj — M k ^ are symmetrical by 
construction. These constant tensors Mf k completely de- 
termine the RG flows. 

The solution for Eq. (@J can be viewed as the trajectory 
of a strongly overdamped particle under the influence of 
the external force Fi in the multi-dimensional coupling 
space. Naively, one might rush to the conclusion that the 
conditions for the existence of a potential requires are, 



dFj dFj 
dgj % 



0, 
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which implies that the tensor Mf k is totally symmetric. 
It is straightforward to check that the RG equations for 
the Q1D systems do not satisfy this criterionS^i^. 

However, under general linear transformations of the 
couplings hi{l) = LijgjQ), the coefficients M? trans- 
form into a new set of coefficients N( , which may be- 
come symmetric. For convenience, we introduce a set 
of matrices [M(k)]ij = M\ to represent the coefficients. 

The symmetric criterion for Nf k — NJj k requires the exis- 
tence of a constant matrix L which satisfies the following 
constraints (for all hi), 



M{k) T = (L T L)M(k){L T L) 



(6) 



where superscript T means transpose. In general, there 
is no guarantee why the strongly over-determined con- 
straints would allow a solution for L. In fact, it is a 
nontrivial task to just prove/disprove whether the de- 
sired linear transformation L exists. Surprisingly, for the 
Q1D systems, the hunt for the solution greatly simpli- 
fies if we formulate the problem in terms of Majorana 
fermions. The desired transformation becomes diagonal, 
Lij = riSij, where r, is a set of rescaling factors and leads 
to the totally symmetric coefficients, 
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So the search for the potential is now nailed down to 
find a set of rescaling factors fj in the Majorana rep- 
resentation. In later section, we demonstrate how to 
construct the RG potential explicitly in the Majorana 
representation. In fact, one can also construct the po- 
tential V(hi) directly from the RG equations for doped 8 
and half-filled 10 ' 16 Q1D systems. Both approaches lead 
the the same result and the detail work will be described 
elsewhere. 

Before we leave this section, it is important to discuss 
a special set of analytic solution of Eq. which is closely 
related to the scaling Ansatz in Eq.^ Suppose the initial 
values of the couplings are <?j(0) = Gig(0), where g(Q) = 



U <C 1 and Gi are order-one constants satisfying the 
non-linear algebraic constraint, 



Gi = Ml k G 3 G k 



(8) 



It is straightforward to show that the ratios between cou- 
plings remain the same and the complicated equations 
reduce to single one, 



dg 
dl 
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For repulsive interaction U > 0, the above equation can 
be solved easily g(l) = l/(ld — I), where the divergent 
length scale Id = 1/U. Note that this implies the ratios 
of different couplings remain fixed in the RG flows, 
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These special analytic solutions are referred as "fixed 
rays" because the ratios of the renormalized couplings 
remain fixed along the flows. One immediately notices 
that these special set of solutions are nothing but the 
peculiar Ansatz found in the numerics. As explained in 
the introduction, if the RG potential exists, these fixed 
rays are the asymptotes of the "valleys/ridges" of the po- 
tential profile and capture the ultimate fate of RG flows 
completely. 



III. QUASI-1D LADDER 

Since the readers might not be familiar with Majorana 
representation, it is worthwhile to write down the Hamil- 
tonian in terms of the familiar field operators for electrons 
first. To construct the RG potential, we first write down 
the general interacting Hamiltonian in weak coupling. To 
be concrete, we would take the iV-leg ladder as the ex- 
ample. However, the general framework developed here 
can be applied to more general Q1D systems. In weak 
coupling, it is natural to work on the band structure first. 
Suppose the chemical potential cuts through Nt flavors 
of bands at different Fermi momenta k Fi with velocities 
Vi. In the low-energy limit, the electron operator can be 
decomposed into chiral fields, 



(ii) 



where a is the spin index, and i stands for the band index. 
The effective Hamiltonian density is simply a collection 
of Nf flavors of Dirac fermions, 



^0 = W Ria (-iVid x )i)) Ria + ipl ic ,(iVid x )tp Li 



(12) 



The summation over the band index i runs through all 
Nf flavors. In general, the velocities are different and can 
not be eliminated by rescaling the space-time coordinates 
(t, x) as in strictly one-dimensional systems. 

Writing down all possible interactions is more laboring. 
By dimensional analysis, the most relevant interactions 
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are the marginal four-fermion vertices. It turns out that 
the allowed vertices can be group together elegantly in 
terms of SU(2) currents. For convenience, we introduce 
the following SU(2) scalar and vector currents, 



jpn = ^i, a ^ a ptpp 3 0, (14) 

Ipu = -ip Pia e a /3ipPi0, (15) 

Ipij = ^IpPid^afilpPifi, (16) 



where ei2 = —£21 = 1 is the anti-symmetric Levi-Civita 
tensor and a are the Pauli matrices. The factor 1/2 
ensures the currents satisfy the conventional SU(2) com- 
mutators. 

At half filling, the particle-hole symmetry pairs up the 
Fermi momenta with the relation 

k Fi + k Ft = 7T, (17) 

where 1 = (N + 1) — i. Other than the above relation, 
different Fermi momenta are generally incommensurate. 
In loose terms, the Fermi surface (points, in fact) is only 
nested between the pairs of bands i and t, with the nest- 
ing condition i + i = N + 1. After some algebra, it is 



straightforward to show that all momentum-conserving 
vertices can be expressed in terms of products of the 
SU(2) currents, 

"WW — rP J J _ rP J .J 

/L int ~ Rij u Lij 'hj" Rii " Li) 

+ fijjRuJhjj ~ fijJrtii ' J Ljj 

"F SijJRijJhji ^ijJ Ri 3 ' ^ L Jt' (1^) 

The couplings £y and /y denote the Cooper and forward 
scattering between bands i and j. The strange vertex 
Sij, involving four different bands, arises from the pair- 
ing up of Fermi momenta, as described in Eq. 1171 The 
superscripts p, a label the charge and spin sectors of the 
couplings respectively. 

Because fa, cu describe the same vertex, to avoid dou- 
ble counting, we choose the diagonal piece of the forward 
scattering amplitude to vanish, i.e. fa = 0. The same 
reason leads to §u — = Su. While it is not obvious at 
this point, the choice of signs for the scalar and vector 
couplings in Eq. ^| is such that they are all positive for 
the repulsive on-site interaction. 

In addition to the vertices which conserve momentum 
exactly, there are also Umklapp interactions which con- 
serve momentum only up to reciprocal lattice vector, 
AP = ±2tt, 



H 



(2) 



In,,) 



IV" 



2 (iij L]S + iy Rii ) - -f(ii t ■ i Lj3 + 1{ • i RU ). 



(19) 



Again, since Uu,Wu describe the same vertex, we set 
u>u = = wa to avoid double counting. Note that both 
kinds of vertices involve four different bands in general. 

Away from half filling, the relation in Eq. El is n0 
longer valid. As a result, the couplings Sij,Uij,uiij all 
become irrelevant because the fast oscillating phase arise 
from the finite momentum associated with the vertex. 
Therefore, we are left with the familiar Cooper and for- 
ward scatteringi^ in the doped A^-leg ladder. 

All vertices discussed in above share a unique feature. 
While the vertices in Eq. ^] and ^] may involve four 
different bands, there are at most two different veloci- 
ties associated with each vertex because Vi = vt, i.e. the 
velocities in each vertex always appear pairwise. This 
seemingly useless feature turns out to be strong enough 
to guarantee the existence of the RG potential when 
the Hamitonian is re-expressed in terms of Majorana 
fermions. 



IV. MAJORANA REPRESENTATION 

Now we switch to the Majorana fermion basis and con- 
struct the RG potential explicitly. Without interactions, 
the band structure in low-energy limits is described by 
Nf flavors of Dirac fermions with different velocities in 
general. Each flavor of Dirac fermions can be decom- 
posed into two Majorana fermions. Combined with spin 
degeneracy, the 4Nf flavors of Majorana fermions are de- 
scribed by the Hamiltonian density 

Ho = VRa(-iVadx)ilRa + ^La (iv a d x )rj La , (20) 

where v a denotes the Fermi velocity for each flavor. 

In general, a single vertex involves four different Fermi 
points, which generally would have four different veloc- 
ities. However, the momentum conservation in weak 
coupling gives rise to the interesting constraint that the 
Fermi velocities must equal pairwise, as demonstrated ex- 
plicitly in previous section. The interacting Hamiltonian 
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FIG. 1: Forward and backward vertices. Bold lines represent 
right-moving Majorana fermions while dashed lines stand for 
the left-moving ones. 



in terms of the Majorana fermions take the form, 

Hint = F(a,a';b,b')7] Ra i] Ra ,ri Lb ri Lbl 

+ B(a,b;a',b')r] Ra r] Rb r] La ,r] Lb ,, (21) 

where summations over allowed indices are implied. We 
emphasize that the allowed interactions might involve 
four different Fermi points labeled by a, a', b, 6', but only 
two different velocities v a — v a > and Vb — vy appear in a 
single vertex. By direct comparison, it is clear that the 
F— vertex include fij and my and the B— vertex cover 
Cij 7 $ij and Uij ■ 

There are also other kinds of interactions allowed by 
momentum conservation, 

H c = R(a,a';b,b')r] Ra r] Ra ,ri Ri r] Rb , 

+ L(a,a' ,b,b')r) La r] La ,ri Lb r) Lb ,. (22) 

The scaling dimensions of these vertices are (A H , A L ) = 
(2,0), (0,2), while the vertices in Eq. 121|) have scaling 
dimensions (1, 1). Since the renormalization comes from 
loop integrations, only vertices with scaling dimensions 
differed by (n,n), where n is an integer, would renor- 
malize each other. As a result, the chiral vertices in 
Eq. (1221 remain marginal and only renormalize the cor- 
responding Fermi velocities. Since the corrections only 
show up at two-loop order, we would ignore their con- 
tribution here. The pairwise-equal Fermi velocities in 
Eq. (|21|l make the classification of all vertices fairly sim- 
ple as shown in Fig. The names come from the fact 
that the forward-type (F-) vertices include the usual for- 
ward scatterings while the backward-type (£>-) vertices 
include the backward scatterings. 

To obtain the flow equations, we need to integrate 
out fluctuations at shorter length scale successively. The 
most convenient approach is by the operator product ex- 
pansions (OPE) of these vertices which form a close alge- 
bra. The detail techniques to compute the renormalized 
interaction STLr can be found in Ref. 

To one-loop order, the renormalization of the bare 
couplings come from four types of diagrams FF — > F, 
FB — ► B, BB — > F and BB -> B, shown in Fig. Let 
us start with the first type of loop diagrams, FF — > F in 
Fig. Ola). The OPE of Majorana fermions can be com- 
puted straightforwardly and the mode elimination leads 







FIG. 2: Four different diagrams to the one-loop order. (A) 
FF -> F (B) FB ~* B (C)BB -» F (D) BB -> B. Notice that, 
only in the fourth diagram, there are three velocities involved 
while only two velocities are involved in all other diagrams. 



to the renormalized Hamiltonian density, 



Sn H = F{a,a";b,b")F(a",a';b",b') 



'■ h" U'^RaVRa'VLbVLb' 



dl 



n{v a + v b ) 

= dF(a,a';b,b')r] Ra r] Ra ,r] Lb r] Lb ,, (23) 

where dl = In 6 is the logarithmic length scale. Here we 
have used the fact that v a = v a ti = v a > and Vf, — Vyi — 
vy ■ The factor l/iv(v a + Vb) arises from the product of 
propagators with opposite chiralities and different veloc- 
ities. 

Introducing a simple rescaling of the original couplings 
according to their associated velocities 



F(a,a';b, b') 



1 



2-K^V a Vb 



F(a,a';b,b'), (24) 



the RG equation is 
dF(a,a';b,b') 



dl 



= lab F( ai a";b,b")F(a",a';b",b') + 



(25) 

where y a b — 2^/v a Vb/ (v +w&). Repeating similar calcula- 
tions, the RG equation for F(a, a"; b, b") contains a term 
r y a bF(a,a';b,a')F(a" ,a';b" ,6) and similar result for the 
renormalization of F(a" ,a';b" ,b'). Therefore, the flow 
equations can be derived from a potential, 

V(F) = - lab F{a, a 1 - b, b')F(a, a"; b, b")F(a", a'; b"b'). 

(26) 

The RG potential for the (FBB) and (BBB) cases 
shown in Fig. Efb-c) and [3d) can be constructed in a 
similar fashion. Finally, combining all contributions to- 
gether, the weak-coupling RG flows for Q1D systems are 
described by the potential 

V(F, B) = -B(a, b; a', b')B(b, c; 6', c')B(c, a; c', a') 
— JabB(a, b'; a, b)B(a", b'; a', b )F(a, a"; b, b") 
- lab F{a, a'- b, b')F(a, a"; b, b")F{a" ', a'; b"b% (27) 
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Summations over all allowed indices are again implied. 
The merits to use Majorana representation enables us to 
construct the RG potential explicitly. Note that if one 
starts from a 'wrong' basis, it is far from trivial to realize 
the fact that the non-linear flows can be derived from a 
single potential. However, in the Majorana representa- 
tion, the linear transformation to the potential basis is 
diagonal Lij = rj&y. After appropriate rescaling of cou- 
plings, the explicit form of the RG potential is derived. 
We have checked for the doped and half-filled Q1D sys- 
tems and found all potentials agree with Eq. ljT7l) . 

There is one loose end about the rescaling factors. For 
most physical systems, the rescaling factor is slightly 
more complicated than (2-K^v a Vb)~ x ■ the subtlety arises 
from the degeneracies of the couplings imposed by phys- 
ical symmetries. This is best illustrated by the following 
simple example. Consider the RG equations for three 
couplings gt, where i — 1,2,3, 

-df-l^—Wk- (28) 

jk 

Since \tijk\ is totally symmetric, the corresponding RG 
potential is V(g) — 5132.93- Suppose the system has some 
symmetry, such as U(l) symmetry for charge conserva- 
tion, and the couplings are degenerate 32 = .93- The RG 
equations are simplified, 

d 9i 2 d 92 , on , 
-^=9 2 , ^f=9i92. (29) 

It is straightforward to show that we need to perform 
a rescaling transformation, (Tii,/^) = (<?i>\/2<?2) to ob- 
tain the potential V(h) — hi(h2) 2 /2. In fact, for cou- 
plings with n-fold degeneracy, an additional rescaling fac- 
tor y/n is necessary to bring them into the potential basis. 
Therefore, the total rescaling factor is 

n = ± x [^, (30) 

27T V V a Vb 

where rij is the degeneracy number of the coupling g il 
with Fermi velocities v a and i>b- 

V. CONCLUSIONS AND DISCUSSIONS 

So far, we have shown the existence of the RG poten- 
tial by explicit construction, which proves the the widely 
used Ansatz in Eq. ^ In fact, the asymptotes of the RG 
flows are governed by the special set of fixed-ray solu- 
tions. This in turns explains the simplicity of the phase 
diagram, even though the RG equations are rather com- 
plicated. 

The absence of exotic fates of the RG flows is also found 
in earlier work on 2D melting theory^ or the Kondo 
related problems^. However, the simplicity of RG flows 
in these systems is not quite the same as described here. 
Since the number of marginal couplings in these problems 



is few, analytic solution of the flows often show that it 
is possible to define some conserved quantity associated 
with each flow lines. This is where the simplicity comes 
from. On the other hand, we have also looked into these 
well-known flows to check whether they can be derived 
from a single potential. It is not too surprising that this 
is indeed the case because the requirement of potential 
flows loose up quite a bit when the number of couplings 
is small. 

With the help of non-perturbative Abelian bosoniza- 
tion, it is not essentially important whether the RG flows 
can be cast into potential form beyond one-loop order. 
However, it remains an interesting and open question at 
this moment. Note that the coefficients of the one- loop 
RG equations are unique, protected by the leading loga- 
rithmic divergences. The next order calculations bring in 
lots of complications and subtleties, including the non- 
universal coefficients in the RG equations, velocity renor- 
malization and so on. It is not clear at this moment 
whether it is even sensible to pursue the RG potential 
beyond one-loop order. 

Another interesting issue concerns the connection 
between the potential V(hi) and the Zamolodchkov's 
c- function C'(gi) of (l+l)-dimensional systems with 
Lorentz and translational symmetries^. A generic Q1D 
system we studied here has neither Lorentz invariance 
(due to different Fermi velocities) nor translation symme- 
try (due to Umklapp processes). While both V and C are 
non-increasing along the RG flows, the exact relation be- 
tween them remains unclear at this point. We emphasize 
that the existence of a non-decreasing function C along 
RG flows only implies that dC/dl = (<9C/%) ■ (dg t /dl) < 
and is not strong enough to show that the flows can be 
derived from a potential. Thus, the potential flows are 
closely related to the c-theorem but they are not equiva- 
lent in general. In addition, we do not know any easy gen- 
eralization of c-theorem that does not rely on the Lorentz 
and translational symmetries. However, one can easily 
check that in the special limiting case where Lorentz and 
translation symmetries are restored, the C-function in- 
deed coincides with the potential we find. This indicates 
that there may be a general form of c-theorem waiting to 
be discovered. 

In conclusion, we have shown that the RG transfor- 
mation for Q1D systems in weak coupling is described 
by potential flows. Therefore, neither chaotic behaviors 
nor exotic limit cycles could occur. The different Fermi 
velocities and the degeneracies imposed by physical sym- 
metries give rise to non-trivial rescaling factors, which 
hinder this beautiful structure behind the RG transfor- 
mation. The explicit form of the potential is obtained 
after appropriate rescaling of the couplings in Majorana 
basis. 
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